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We report the plasmon dispersion characteristics of intrinsic and extrinsic armchair graphene 
nanoribbons of atomic width = 5 using a pz-oih\idX tight binding model with third-nearest- 
neighbor (3nn) couphng. The couphng parameters are obtained by fitting the 3nn dispersions to 
that of an extended Hiickel theory. The resultant massive Dirac Fermion system has a band gap 
Eg ^ 64 meV . The extrinsic plasmon dispersion relation is found to approach a common dispersion 
curve as the chemical potential /i increases, whereas the intrinsic plasmon dispersion relation is 
found to have both energy and momentum thresholds. We also report an analytical model for the 
extrinsic plasmon group velocity in the ^ ^ limit. 



I. INTRODUCTION 

Graphene exhibits massless Dirac Fermions^^ with semi- 
metallic behavior, for which the collective carrier modes 
in the form of plasmons have been a topic of stud}^ — . 
When graphene is patterned in the form of an armchair 
graphene nanoribbon (acGNR)ii — , two important devi- 
ations occur. First, the acGNRs develop a band gap and 
second the dispersions do not remain linear anymore and 
hence the electrons and holes behave as massive Dirac 
Fermions. 

While acGNRs of atomic widths TV of mod{N, 3) = 0, 1 
exhibit significant band gap opening irrespective of the- 
oretical model, acGNRs with mod{N,3) = — 1 have zero 
band gap and massless dispersion within the continuum 
and the first nearest-neighbor p^-orbital tight binding 
(Inn pzTB) mode\^^^^. One has to use more de- 
tailed methods like density functional theory (DFT)^ 
extended Hiickel theory (EWT)^^ or beyond Inn TB 
model to get a more detailed band structure. Although 
there are quantitative differences amongst these methods, 
nonetheless qualitatively these methods converge upon 
massive Dirac Fermions with a band gap opening for 
mod{N,3) = — 1 acGNRs. The band gaps predicted by 
EHT for these acGNRs are of the order of a few tens of 
meV for extremely narrow ribbons, and decreases as the 
width of the nanoribbon increased. 




FIG. 1. (Color online) Ball and stick model of the hydro- 
genated acGNR5. x- (y-) is the longitudinal (transverse) di- 
rection. Atomic visualization is done by Hiickel- NV—. 

In this paper, we examine the plasmon dispersion 
in intrinsic and extrinsic acGNR with atomic width 



N = 5 (acGNR5) by using a third-nearest-neighbor (3nn) 
pzTB model, benchmarked with EHT, within the random 
phase approximation (RPA) as discussed in Sec. II. We 
discuss the plasmon dispersion results in Sec. HI followed 
by the conclusions. 



II. THEORETICAL MODEL 

The unit cell for a hydrogen passivated acGNR5 is 
highlighted in Fig. [T]that contains 10 carbon and 4 hy- 
drogen atoms. The unit vector is given as a = = 
SaccX^ where acc = 1.42 A is the carbon bond length. 
The pzTB Hamiltonian of the unit cell is a 10 x 10 ma- 
trix containing 3nn couplings. We transform the real- 
space Hamiltonian to the reciprocal space H{k) to cal- 
culate the eigenvalues Ei{k) and eigenfunctions cf^\k) 
for the eigenstate i = 1, 2, . . . , 10, where i is the band 
index and a represents the atomic location. The band 
index ranges from i = 1 {i = 1) corresponding to the 
lowest-lying conduction (highest-lying valence) band and 
i = 5 {i = 5) corresponding to the highest-lying con- 
duction (lowest-lying valence) band. One finds that the 
electron-hole symmetry is broken due to finite 2nn and 
3nn couplings. 

The band structure for acGNR5 is shown in Fig. [2l 
3nn tight-binding parameters (£^07^07^17^2) for the fit of 
the acGNR5 nanoribbon to the EHT— are reported in 
Table HI Parameters are obtained by fitting the top three 
valence bands and bottom three conduction bands of the 
EHT data to the 3nn pzTB band structure at 51 k-points 
uniformly spaced across the Brillioun zone. Fitting is 
accomplished using a least-squares algorithm and no ge- 
ometric relaxation of the bond lengths is incorporated. 
This set of hopping parameters agrees well with Ref 
for the to and ti parameters. However, the 3nn hopping 
parameter (^2) we report is significantly smaller, due to 
the smaller gap predicted by EHT^^-^^ when compared 
with DFT results^. 

The band structure computed using the 3nn Hamil- 
tonian with hopping parameters from Table |I] is shown 
in Fig. [21 The carrier dispersion characteristics of the 
z = 1, 1 (valence, conduction) bands in this system show 
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FIG. 2. Carrier dispersion for acGNR5 calculated using the 
3nn tight binding Hamiltonian with hopping parameters given 
in Table [H Panel (a) shows the complete 10-band structure, 
and panels (b) and (c) show progressively more detail in the 
dispersion near the gap of Eg ^ 64 meV over two different 
ranges of momentum. In addition to the gap, asymmetry be- 
tween the conduction and valence bands has been introduced 
by the non-zero 2nn and 3nn coupling parameters. 



TABLE I. 3nn tight Binding Parameters for the acGNR5 
nanoribbon obtained using a fit to EHT^^ data. 
Eq 0.11031 eV onsite energy 
to -2.69341 eV Inn hopping parameter 
ti 0.02201 eV 2nn hopping parameter 
t2 -0.03225 eV 3nn hopping parameter 



a finite gap of Eg ^ 64 meV and represent a massive 
Dirac Fermion system with a dispersion characterized by 
the relation: 
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v&^^^)^TS^^^ (1) 



where vpi is the Fermi velocity for the z = 1, 1 (va- 
lence, conduction) bands, and the + (— ) sign is chosen 
for the conduction (valence) band. The band gap Eg cor- 
responds to a relativistic rest mass of the massive Dirac 
Fermion system of mo = Eg/2vp-. 

To compute the plasmon dispersion in the random 
phase approximation (RPA) for the nanoribbon with 3nn 
Hamiltonian, we follow the procedure outlined in Ref.^^. 
Due to large energy differences and small electronic wave- 
function overlap integrals at g ^ for both Inn pzTB in 
Ref and 3nn pzTB in this paper, we use a two-band 
dielectric function including only the i = 1 conduction 
and i = 1 valence bands to study the plasmon dispersion 
relation. However, some of the details of the 3nn pzTB 
model are different, which we discuss next. 

In the RPA expression for the interband polarizability, 
electronic wavefunction overlap integrals between states 
in the two bands at momenta k and k' = k -\- where q 
is the plasmon momentum, play a significant role. The 
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FIG. 3. 3nn tight-binding overlap integral computed for 
(a) the interband transition between conduction and valence 
bands, and for (b) the intraband transition of the conduction 
band along the line k' — k -\- q. In each panel, a family of 
12 curves is shown for < < qmax where qmax = llAg^, 
in steps of where = 7v/{800d). In (a) the overlap for 
^ = is identically 0, and the overlap for q = qmax has a 
maximum of approximately 0.66, whereas in (b) the overlap 
for ^ = is identically 1, and the overlap for q — qmax has a 
minimum of approximately 0.76. The vertical bars show the 
bounds —qmax < k < 0, corresponding to the region where 
1(1; /c^ = k -\- qmax\^] k) \ — 1 in the continuum model. Outside 
of this range, \{l;k' = k -\- qmax\'^'i k)\ = for the continuum 
model. The conduction-conduction band has the opposite 
symmetry in the continuum model. 



polarizability is written as 

Umniq,^) 
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\{n',k'\m',k)\^ 



(2) 



where m and n are band indices, = 2 is the spin de- 
generacy, Lx is the sample length, k is the momentum 
of the initial state, k' = k ^ q \^ the momentum of the 
final state, and f{E) = 1/[1 + e^^-^^/^^^] is the Fermi- 
Dirac distribution function with chemical potential fi and 
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Boltzmann constant ks, where T is the temperature in 
K. h is the reduced Planck's constant and 7^ is a smah 
number. We consider intrinsic acGNRs with the chemical 
potential /i = 0. In Fig. [3l we illustrate several of these 
overlap integrals as functions of k. It should be noted 
that the the overlap integral is no longer confined in the 
region defined by sign(/c/c') = —1, as it is in the contin- 
uum model. Rather, the overlap integral is non-zero well 
beyond this range. However, direct transitions at g = 
are still forbidden. The broadening of the overlap inte- 
gral beyond the hard boundaries of the continuum model 
indicates that there is a coupling to free-carrier states for 
collective modes at all nonzero and as a result, plas- 
mons in this system are Landau damped. 

Following Refjiii, we calculate the longitudinal dielec- 
tric function for acGNRs in the RPA. In the RPA, the 
dielectric matrix for acGNRs can be written as^^: 




FIG. 4. Extrinsic plasmon dispersion relations for acGNRS. 
Dispersion curves are calculated for a range of chemical po- 
tentials /i. 



{q)Umn{q,(^) (3) 

where VijmniQ) is the Coulomb matrix element in one 
dimension, UmniQ^ ^) is the polarizability of the acGNRs, 
and z, j, m, and n are the band indices. Non- trivial 
solutions to the field equations require: 



det [eijmn{q,(^)] = 



(4) 



Intrinsic Plasmons: In the two-band approximation 
for intrinsic acGNRs at T = 0, the self-polarizabilities of 
the i = 1, 1 bands are given as, Iiji{q^uj) = Ilii{q^uj) = 
0. Further, symmetries in the acGNRs require^^i^ that 
the Coulomb matrix elements '^iiii(^) = y i i t(^) = 
'^1 T 1 t(^) ~ '^iTTi(^)- This result gives the dispersion 
relation of the collective (plasmon) state in the 2-band 
approximation by simplifying Eq. 3] as follows: 



(5) 



We compute the Coulomb matrix elements vj ^ (q) as 
described in Ref.= using the p^-orbital wavefunction lo- 
calization parameter w = 1A2^. Solving Eq. Ogives the 
dispersion relation for the collective modes (plasmons) in 
the acGNR. 

Extrinsic Plasmons: The dispersion relation for plas- 
mons in extrinsic acGNR can also be obtained from Eq. 
[4] in the two-band approximation. For a chemical poten- 
tial /i in the i = 1 conduction band at T = 0, states with 
momenta between —kf < k < kf where Ek^i = /i are 
filled, and states outside of this range are empty. For the 
extrinsic case Ilii{q^uj) is no longer 0, and we write the 
plasmon dispersion relation as: 



(1 -^T,l,T,i(^) i^TiiQ^^) ^^iiiQ^^)]) 
X (1 - vi,i,i,i((7)nii((7,a;)) = 



(6) 



Plasmons for negative chemical potentials /i will exhibit 
similar behavior. 
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FIG. 5. Group velocity of extrinsic plasmons as a function 
of the chemical potential fi in the ^ ^ limit. The solid 
points are calculated from the dispersion relation data pre- 
sented in Fig. m The dashed curve is calculated using the 
analytic model for limg^o Vg (/i) discussed in the text with 
^^i,i,i,i(0)/(2eVeo) = 11.1294 and vfi = 8.33219 x 10^ m/s. 



III. DISCUSSION OF RESULTS 

Intrinsic Plasmons: The intrinsic plasmon obtained 
using our formalism exhibits an onset threshold in both 
the q and E dimensions. The q threshold can be under- 
stood with the data presented in Fig. [3l^a). For small val- 
ues of the overlap integral is nearly zero. Because the 
polarizabilities Ilj-^{q^ uj) and n-,^Y are proportional to this 
overlap, the dielectric function never crosses 0, and so no 
collective mode exists. As the overlap gets larger, the di- 
electric function eventually crosses zero and the intrinsic 
plasmon dispersion exists. The threshold in £^ is a result 
of the fact that the polarizabilities are not large enough 
to cause a zero-crossing for small E. As the plasmon en- 
ergy increases above the bottom of the conduction band, 
the resonant enhancement in the polarizabilities causes 
a zero-crossing. Because we are interested in plasmons 
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in the ^ ^ limit, we do not consider the intrinsic case 
further. 

Extrinsic Plasmons: Dispersion relations for plasmons 
in extrinsic acGNR5 computed using the tight-binding 
formalism described above are also plotted in Fig. |4]for 
several values of the chemical potential fi > Eg/ 2 corre- 
sponding to a geometric distribution oi kp- From these 
results, it can be readily observed that the dispersion 
curves have a character for values of the 

chemical potential within a few meV of the band edge 
{ji > Eg/2). Further, as the chemical potential increases, 
the dispersion relation is observed to asymptotically ap- 
proach a limit which corresponds to the plasmon disper- 
sion in a massless Dirac fermion system. 

It is interesting to analyze the behavior of the extrinsic 
plasmon group velocity in the g ^ limit as a function 
of the chemical potential fi. In this limit, the interband 
polarizabilities II-,^y((7, cj) = nYx(<7,cj) = because the 
interband overlap integral (1; /c) = (see Fig. [3fa)). 
Further, the intraband overlap integral (1; k\l] k) = 1 m. 
this limit (see Fig. [3fb)). As a result, the intraband 
polarizability becomes: 

nn(<Z,-) = ;^^^ (7) 

k ^ ' 

where l\E = £^/c+g,i —E^^i. Thus, as ^ 0, the dielectric 
function becomes: 

i-^i,i,i,i(g)nii(g,cj) = (8) 

Solving Eq. [8] for cj, the plasmon group velocity in the 
q ^ limit can then be written: 



Vgikp) = 




where the integral is taken over the filled states be- 
tween —kp and kp- Substituting the relationship be- 



t ween chemical p otential and Fermi wavenumber kp = 
Y^/i^ — (moVp-^)'^ /hvpi into the analytic result for Eq. [9] 
gives the group velocity as a function of chemical poten- 
tial Vg{jii)^ which is shown as a dashed curve in Fig. [5l 



IV. CONCLUSIONS 

In summary, we have computed the plasmon disper- 
sion for an acGNR5 nanoribbon using a 3nn tight-binding 
model. This nanoribbon represents a massive Dirac 
Fermion system. Hopping parameters for the model were 
obtained by fitting the 3nn band structure to band data 
obtained from an EHT calculation. The intrinsic plas- 
mon dispersion relation obtained exhibits a threshold in 
both q and E. The extrinsic plasmon dispersion rela- 
tion obtained follows the qy/V{q) dependence expected 
in ID systems for values of the chemical potential near 
the band edge (/i > Eg/2)^ and the dispersion relation 
asymptotically approaches one corresponding to a mass- 
less Dirac fermion system as the chemical potential ji in- 
creases. Good agreement between the group velocity of 
these plasmons in the q ^ {) limit and an analytic model 
based on the behavior of the polarizabilities as g ^ is 
obtained. Finally, we note that some damping of these 
plasmons may be expected to occur from plasmon scat- 
tering to free electron states due to the nature of the 
relevant overlap integrals. 
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